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Static solutions of the electro-gravitational field equations exhibiting a functional relationship 
between the electric and gravitational potentials are studied. General results for these metrics 
are presented which extend previous work of Majumdar. In particular, it is shown that for any 
solution of the field equations exhibiting such a Weyl-type relationship, there exists a relationship 
between the matter density, the electric field density and the charge density. It is also found that the 
Majumdar condition can hold for a bounded perfect fiuid only if the matter pressure vanishes (that 
is, charged dust). By restricting to spherically symmetric distributions of charged matter, a number 
of exact solutions are presented in closed form which generalise the Schwarzschild interior solution. 
Some of these solutions exhibit functional relations between the electric and gravitational potentials 
different to the quadratic one of Weyl. All the non-dust solutions are well-behaved and, by matching 
them to the Reissner-Nordstrom solution, all of the constants of integration are identified in terms 
of the total mass, total charge and radius of the source. This is done in detail for a number of 
specific examples. These are also shown to satisfy the weak and strong energy conditions and many 
other regularity and energy conditions that may be required of any physically reasonable matter 
distribution. 
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I. INTRODUCTION 

The static field has played a central role in our understanding of Einstein's theory of gravitation. Not only does this 
area lie closest to Newtonian gravitation and classical potential theory, but the existence of a timelike hypersurface 
orthogonal Killing vector brings about many simplifications in the field equations which allow one to integrate them 
exactly in numerous situations. 

In particular, the (744 component of the metric tensor plays the part of a gravitational potential. In the presence 
of an electric field a similar electric potential (f) exists, and it was by considering fields in which there is a functional 
relationship between the two potentials that Weyl Q discovered a class of axially symmetric solutions. We will refer 
to any static solution of the field equations with such a functional relationship as being of Weyl-type. 

Weyl considered only pure electric fields (that is, fields in which the electric field is the only source of the energy 
momentum tensor). For these fields Weyl found that, in the presence of axial symmetry, the functional relationship, 
if it exists, must be of the form 

g44 = A + B(j) + ^\ (1) 

Later work by Majumdar g showed that, even in the absence of axial symmetry, this quadratic restriction persists. 
Majumdar also discovered that, if this quadratic relationship is a perfect square 

344 = (C + <j>f , (2) 

then the spatial sections are conformally flat and one can interpret the space-time as describing the exterior field of 
a static ensemble, where the gravitational attraction and electric repulsion exactly balance. These pure electric fields 
were independently discovered by Papapetrou y] . We will refer to any static solutions of the field equations with a 
perfect square functional relationship (0) as being of Majumdar-type. 

The case of Weyl-type fields in which there are sources in addition to an electric field has received some attention. 
Majumdar discovered solutions which represent the interior of a distribution of charged dust, held in balance by 
gravitational and electric forces. Here again the functional relationship is a perfect square (y) and the matter pressure 
vanishes, while the matter density is numerically equal to the charge density. In addition, the total mass and total 
charge of the space-time are found to be numerically equal. 

Other than work done on charged dust, where the matter pressure vanishes (for example Das Q, Giirses [g |g] 
Raychaudhuri 0|) little has been done to investigate Weyl-type fields in the presence of matter. The question we 
address in this paper is whether the functional relation in the presence of matter sources (in particular, matter 
pressure) must be Weyl's quadratic one (M. The answer to this is no, and we exhibit physically reasonable solutions 
to the field equations where the functional relationship can be in a range of forms. 

This paper is organized as follows: in the next section we give the background, and prove a number of general 
results for Weyl-type fields, with or without spatial symmetry. We show, generalising Majumdar's results for charged 
dust, that for Majumdar-type fields, the "active gravitational matter density" is numerically equal to the charge 
density. We show that for asymptotically fiat distributions of charged matter of Majumdar-type, the total mass and 
total charge, as measured at spatial infinity, are numerically equal. Finally, we show that any bounded non-singular 
perfect fluid which is of Majumdar-type consists of dust (i.e. is pressure free). 

Sections 3 to 5 looks at the spherically symmetric case. We present explicit Weyl-type generalisations of the interior 
Schwarzschild metric, in which the functional relation takes various forms. As far as the author is aware, these are 
the flrst explicit solutions in which the functional relationship deviates from Weyl's original quadratic one. In section 
4 a number of these are then investigated as regards their matching conditions with the exterior Reissner-Nordstrom 
solution. This allows us to identify the constants of integration in terms of the total mass, charge and radius of the 
sphere. In section 5 it is also found that many satisfy appropriate regularity and energy conditions. 

II. BACKGROUND AND GENERAL RESULTS 

A space-time (A/, g^^) is static if there exists a timelike hypersurface orthogonal Killing vector ^^ on M. Here, and 
throughout, Greek letters will take values 1 to 4. With a suitable choice of co-ordinates (a;*,t) the metric takes the 
form 

ds"^ = gijdx'-dx^ - e'^dt^, 

where both the spatial metric g^j and e" only depend on a;*, for i,j=l to 3. 



Einstein's field equations say that the Ricci tensor R^i, of g^i, is related to the energy-momentum tensor T^^, by 

where R is the scalar curvature R — g^'^ R^^. In the case of charged matter, the energy momentum tensor is made 
up of two parts, the matter tensor Af^i^, and the electromagnetic energy tensor E^j^^: 

The electromagnetic energy tensor is determined by the (skew-symmetric) electromagnetic tensor F^^, 

ifi-R — F "F — -F aF^l^n 
which satisfies Maxwell's equations 

V^F^"" = J", 

V[„F^,] = 0. (3) 

Here V is covariant differentiation with respect to the Levi-Civita connection of g^^. J'^ is the 4-current density 
vector: 

where V^ is the 4-velocity and a the proper charge of the charged matter. Equation (0) is the integrability condition 
for the existence of the local 4-potential 0^ : 

F^u = V^ (/),. - V,. (t>fj,. 

We will consider only pure electric fields 0^ = (0,0,0,0) (the magnetic field is assumed to vanish), noting, however, 
that a theorem of Bonnor |^ allows one to generate pure magnetic fields from pure electric fields. 
For static pure electric fields, the non-vanishing components of the electromagnetic energy tensor are 

KEij = -2e""(/),j(?!),j -(-e""(/)'''0,fc5jj, 

A comma here indicates differentiation, and (j) is the electric potential: 

Fa = (f>,i Fij ~ 0. 

Note that cj) is only defined up to an additive constant (/) ^ (p + C. 
In addition, the 44 component of the Ricci tensor is determined by u 

i?44 = -ie"(Aco + i|dwn, 

where A is the covariant laplacian Aw = g''^LU\ij, a stroke representing covariant differentiation with respect to the 
Levi-Civita connection of the spatial metric gij . 
Maxwell's equations reduce to 

A</.= iw'V„+iKae5". (4) 

The Einstein equations are 

i?y - 2e-^cj),^ cj),j -nM,j + {\kM^ - e-'^\ddp\'')g,, (5) 

-ie"(Atj + i|da;|2) = -\n{A4 - Mt)e^ - |d0|2 (6) 

R = kM^. (7) 

Equations (^), (||), (||) and (^ are the equations governing the static electrogravitational field. A Weyl-type field is 
an electric field satisfying the above equations where 344 and 4> are functionally related: uj = uj{4>). 



Theorem II. 1 For any Weyl-type field the following relationship holds between the gravitational matter density, the 
electric field density and the charge density: 



Mt - Mt 



d^e- 



Et ^ i^.. 



(8) 



Proof: 

The result comes from comparing the 44 field equation (|^) with Maxwell's equation (p. Assuming that lo = u}{(i)) 
we have that 



duj 



\M' = ( S 



^'''4',k- 



duj 



M' 



Aw 



duj 



Ao 



d^LO 



M' 



Substituting these in (g) we get 



2'^ 



duj . 

^ A<; 



and then using Maxwell's equation 



1^1^ 



dw 



^"ijJ , , ,,n . I duJ 

372 \d^\ +2\-n 



duj 



372+37 m 



^_i«(Affc-M4V- 



= -\n{A4 - Mt)e- - 



Rearranging this equation gives the result. □ 

This theorem generalises Majumdar's result g]: 
Corollary II. 2 For a pure electric field (M^^, = and a = 0), the only possible form of a functional relationship 

LU = Lij{(j)) is 



e'^ = A + B(f) 



A, B constant. 



We also have the following two corollaries: 
Corollary II. 3 For any Weyl-type field satisfying the relation e^ = A + B<j) 



Mt - Mi 



1, 



where A = (-j) 



A. 



Corollary II. 4 For a Majumdar-type field (i.e. one where e" = (-j ± (j))^ ) the following holds: 



Mt - Ml = ±(T. 



Thus the active gravitational matter density and the charge density of a Majumdar-type field are numerically equal. 
A further feature of this balancing in Majumdar-type fields is that the total mass and total charge are equal: 
The total mass of a static asymptotically flat space-time {M,g^^) is given by (Wald Q) 



e 



lim — , 



-V^C.iV^ dA, 



where Sr is a sphere of radius R in the hypersurface orthogonal to the KiUing vector ^'^, A'"'^ is the outward pointing 
normal to Sr which is orthogonal to ^^, and integration is with respect to the volume element induced on Sr by gij. 
On the other hand the total charge is given by (Wald jO) Synge pc|]) 

q= lim -^ / Ff^.N^e dA, 

where f^ is the unit normal to the volume enclosed by 5*^. 

Theorem II. 5 For an asymptotically flat Majumdar-type field m = \q\. 

Proof: 

The choice of co-ordinates we have made implies that the Killing vector is given by 



so that 



and the total mass is given by 



^^ - S^ Cm = -e-54^ 



e^V^ C. = 5^"^- = 5^" S <^-. . 



m= lim -^ / e5" ^ 0,, Ar'^dA 
fl^oo Stt Js^ d(p 



On the other hand the total charge is 

q^ lim ^ [ e~^'^(l),kN''dA. 
fl^oo An Js^ 

For a source satisfying Weyl's condition 

e'^ = A + B(j) + (l3^, 

for which the surface of the body is an equipotential surface, Majumdar yj has shown that this relationship must 
continue to hold outside the source (see also Gatreau and Hoffman [ll|). In any case, we will assume that Weyl's 
condition holds throughout the space-time. Then 



m= hm — / ^ (p,k N^dA 

fl^oo Stt Js„ y'A + B(t) + 02 

= lim -!- / {B + 2<j))e-i'^(j),kN''dA. 



'Se 

We assume that the space-time is asymptotically flat, so that — O(-) and A = 1. Then 

Ft f 

m= lim -- / e^i'^(j),kN''dA 
B.'^oo Stt Jg^ 

^Bq 

2 ' 

For asymptotically flat Majumdar-type fields B = ±2 and so to = \q\ as claimed. □ 

Finally we look at the compatibility of Weyl's condition with a perfect fluid matter source. Thus, we assume that 
the matter tensor is of the form 

Mi^ ^ (p + p)y^K+P<5^, 

where V^ is the 4- velocity of the fluid and p, p are the matter pressure and density, respectively. We will show that 
for a static field, Majumdar's relation cannot hold for any physically realistic perfect fluid unless the matter pressure 
vanishes. This is one of the consequences of the following: 



Theorem II. 6 For any static electrogravitational field in which the matter tensor is a perfect fluid, the following hold: 

(i) If any two of the surfaces of constant 544, (f> or p coincide, then the third also coincides. 

(ii) If the fluid is of Weyl-type, then either the pressure gradient vanishes at the surface of the fluid or the surface is 
an equipotential surface (i.e. (j) and e'^ are constant on the surface). 

(Hi) If the fluid satisfies Majumdar's relation ^, then the pressure satisfies 

p — Ce^ C constant., 

and the spatial slice is conformal to a space of constant curvature. 

Proof: 

(i): To prove this we consider the conservation equation 

V^T,^ = 0. 
For static electrogravitational fields the conservation equation v — i reads 

p,r+^ip + p)uj,i = (Je-^'^(l),i. (9) 

From this we see that if any two of p,i, uj,i and (/),.i are coincident, then all three are. This proves (i). 

(ii): We prove this by using the Synge-O'Brien junction conditions ||l^. These state that, if E is the spatial surface 
formed by the boundary of the body with normal rji , then 

Gijri^vi-' is continuous across E. 
Now, since the electromagnetic field is continuous across E (see e.g. Ehlers [Esi) 

Mijrfrj-' is continuous across E, 

and hence, for a perfect fluid we flnd that 

p(E) = 0, (10) 

i.e. the surface of the body is a surface of constant (zero) pressure. 

Therefore if 544 and are functionally related and the pressure gradient does not vanish on E, then, by part (i), E 
is an equipotential surface, as required. 



II.4 



tells us that 



(iii): Finally, note that, for a perfect fluid satisfying Majumdar's condition. Corollary 

3p + p = ±cr. 
Plugging this and Majumdar's relation into the conservation equation (B), we flnd that 

Pri = P^n , 

which integrates up to 

P = Ce'^, (11) 

for some constant C. Thus we flnd that 

M,, = Ce'^g,j. 

Now, inserting this in the fleld equations (pi) and making a conformal change of metric gij — e^'^^j , we flnd that the 
Ricci tensor of 5^ is 

Rij = 2C5,y, 
and this is the deflning equation for a 3-dimensional space of constant curvature. □ 



From this we can draw the conclusion that any bounded non-singular perfect fluid satisfying Majumdar's relation 
must be dust (i.e. pressure- free) . This follows from the fact that if C ^ in ( [ll|) then the boundary condition (|l(]| ) 
will mean that at the surface S we must have 

e"(^) = 0, 

which is impossible, unless the space-time is singular at the boundary. Thus the matter pressure must vanish through- 
out. In the next section we look at a number of spherically symmetric perfect fluid solutions with various functional 
relationships between e" and 4>. There we will see an example of this phenomenon of the incompatibility of Majumdar's 
relation with a bounded perfect fluid, unless the matter pressure vanishes. 

III. SPHERICALLY SYMMETRIC SOLUTIONS 

In this section we illustrate how Weyl's assumption 

544 = 544 (</»), 

can be used to obtain exact solutions of the Einstein-Maxwell field equations in the case where neither the matter 
density nor the matter pressure vanish. In order to do this, we restrict our attention to spherical symmetry and 
perfect fluidity (i.e. isotropic matter pressure). 

The case of static spherically symmetric perfect fluids has been investigated by numerous authors and many exact 
solutions have been extracted (see for example Buchdahl |I3], Buchdahl and Land p5|, Klein [|l6|, Tolman [ p7[ , 
Whittaker ||l^). Perhaps the best known is the interior solution of Schwarzschild ||l^, in which the matter density 
is taken to be constant. This neutral solution has been generalised to include electromagnetic charge by Bohra and 
Mehra ||§, Cooperstock and De La Cruz ||l|, Florides Q, Mehra H ^ and Patel and Tikekar ||||. We will also 
present a number of charged versions of this familiar solution in what follows, one class of which contains the solution 
of Cooperstock and De La Cruz as a special case. However, the motivation for considering such solutions here is not 
just to add to this list, but to show that Weyl's condition can lead to exact solutions when the field is neither pure 
electric nor charged dust. Moreover, we will show that the functional relation is not constrained to be Weyl's original 
quadratic one (|l]). 

Our first task will be to formulate the static Einstein-Maxwell field equations when the spatial metric exhibits 
spherical symmetry. The metric for such a distribution can be written in local co-ordinates (r, 9, ip, t) as 

ds^ = e^^'^^dr^ -t- r'^ide'^ + sin^Odip^) ~ e'^^^'^dt^. 

From spherical symmetry and the Maxwell equations we find that the only non-vanishing components of Tj^ are 
(see Florides [p6| ) 

nil = np- —J- KJ4 = -up ^ 



j.t 



kT^ = nTi = kM^ + ^, 



2 



where p and p are the radial matter pressure and density, and Q is the total charge inside a radius r, i.e. 

p{r) = Ml 
p{r) = -Mt 

Q{r) = f /" a{ryei°'dr 
Jo 

Now, the field equations 

when we insert the spherical metric form and the energy momentum tensor, read 



Q 1 e~" 



kM^ 



Q' 



e -" ( \u" - \olJ + hJ-" 



2r 



where a prime denotes differentiation with respect to. r. 

The above equations constitute the most general system of equations for any static charged sphere in general 
relativity. However, we still have six unknowns (p, p, M|, Q, a, to) and only three equations. To make the system 
determinate we will first assume that the fluid is perfect [Ml = M| ~ Mg). In this case the field equations become 



np J ^ ~— ^ —il + ruj ) 



KP+^ 



e -" f io;" - ia'c.' + W^ + 



2r 






^ 1 e~" 
Kp+—- = — -(1 - ra'). 



Following Tolman 



we note that subtracting ( ITSI) from (|l^ ) we get 



die 
dr 



dr 



2r 






(12) 
(13) 
(14) 

(15) 



Our plan of attack for solving the field equations (O) to (11^) will be to introduce two conditions which allow us 
to integrate ( p^ ) fully for u: and a as explicit functions of r. Then we substitute these in ( p^ to (|lj) and determine 
p(r), Q[r) and p{r). 

The first condition we impose is the Schwarzschild condition 



Equation (Rj) implies that 



T'4 = Kp 



Q^ 



B? 



constant. 



(16) 



Plugging this in (nSl) we find that 



-'-R^- 



d_ 

dr 



<^'e "\ ^^-i<^+^)d 



1r 



dr V 2r 



40^ 



(17) 



or multiplying across by 



2w'e'^ 



d 


w'2e"-"- 


sg 


^cj'e'^ 


dr 


^2 


^6 ' 


d 


/w'e^("-")\' 


4Q2e" 


dr 


V ^ / 


r^t^' 








4(<^')' 



(e'^)' ■ 



Now we introduce Weyl's condition uj = lu{(J)), or, more conveniently. 



(e'^). Then we have that 



dr 



In 



/ i(w-a)~ 



LJ 6 2 



or integrating 



In 



jlQ^(u-a) 



= 4 



for some positive constant D. Now, define 

/(e") = exp -4 
Then a further integration yields 



= 4 



d0 



(1 



d(e" 



^(6"^ 



dr 



^.1 

_R2P 



= F-DR^[l- — 



(e-y 



In 2D, 



(18) 



(19) 



(20) 



where F is another constant of integration. 

Our strategy to solve the equations will be to choose a functional relation (p = 0(e'^), integrate ( [l^ to get /(e"^) 
and then integrate (EQ) and invert to find e'^ as a function of r. We will investigate two types of functional relation 
and their corresponding solutions, which we refer to as Class /and Class II solutions. 



A. Class I Solutions 



First we consider fields in which 



Ae" 



B + iC + (t>f 



(21) 



where A, B and C are constants. For A = 1 this is precisely Weyl's original condition (m), and if, in addition, B = 0, 
this reduces to Majumdar's condition (g). Indeed the appearance of A in this equation may at first seem superfluous, 
at least for ^4 > 0, as we have the freedom to rescale our time co-ordinate 



t' = A^t, 



so that 



Ae" 



However, if we initiate this co-ordinate change, we must also change the electromagnetic tensor Ffj_^, and it is easy to 
check that this change of co-ordinates induces a rescaling of 



A^ 

Hence, a change of co-ordinates of this type does not alter the form of the functional relationship, since it merely 
induces a rescaling of the constants B and C. We shall see that in order to integrate the field equation with B y^ 0, 
we will allow A to only take on the values 1, |, § and that each of these leads to significantly different solutions. For 
A j^ 1 we are considering for the first time (as far as the author is aware) a functional relation different from the 
original one introduced by Weyl. 

Taking this relation then and computing /(e") from dl9) we have that 



/(e-) = (Ae- - BY 



and then ( pOJ ) tells us that 



{Ae^ - B) 



F - DR"- 1 






For A = Q this gives the standard interior Schwarzschild solution. For B = this can be integrated up to give e"*^'") 
for any value of A (although we require that A > for the functional relationship to be real). For B ^ we can 
integrate this completely if 



A = ±n,±(n- 



n e N. 



However, in general, we then get an implicit equation for e" in terms of r. The only values of A that lead to an 
explicit expression for e" for B ^ are 1, i, |. 

Thus the logic of the field equations forces us to consider the cases B — { Class la) and B ^ { Class lb) separately. 
We now list the metric components, electric charge, density and pressure of each solution. Recall that for each one 



For brevity, we have defined 



= 1 - 



i?2- 



Class la 



Class lb 



^ - 2- 



$(r) = F - DR^ 1 



r 



2 \ 2 



A^h 



Qir) 
Kp{r) 



e'^ = [(1 - 2A)A^$] ^- 
ADr^ 



±- 



"(1-2A)$ 
3 AD^r^ 



i?2 (1 - 2A)2$2 



1 AD2r2 

'^Pir) = -J^+ (1 _ 2A)2$2 



2^(1-1^)^ 
(1-2A)$ 



A=i 



Q{r) = ± 



= \/2$(r) 



, , 3 7?2^2 



^pir) = -^ 



4 

^2^2 



V2i.fl-- 



Q{r) - ± 



2B cosh^ 


T2^i^)_ 


for 


B>0 


-2Ssinh2 


"tj'^w' 


for 


B <0 


Dr^ 
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Kp{r) 



Kp{r) 



3 D^r^ 



i?2 



D'r' 



i?2J • 



V2D tanh 
V2D coth 



V2 



/or B >0 
for B <0 



A=l: 



Kp{r) 



B coth^ 



B$(r) 



Kp{r) 



for B>0 

-Btan2 [V^B$(r)] /or B <0 
±^/BDr^cosech[^/B^] for B>0 
±^/^BDr^scc[^/^B(^] for B<Q 

■^- BD'^r'^cosech'^[y/B<5] for B>0 
-4 + BD'^r'^ sec2 [V^B $] /or S < 



-^ + BD'^r'^cosech'^lVB $] 
- -^ - BD2^2 gg(,2 j^/z;g ^] 



W 



for B>0 
for B <0 



A^i 



2B^ 



F - DR^ 1 



352 



F- DR^ (1- if)' 



, 3DBr^ 
QW=±3S2<j,2_2 

3 9B^D^r^ 



i?2 (3B2$2 _ 2)2 



1 9B^D^r^ 



41? 1- 



i?2 (3B2$2_2)2 (3B2$2_2)$' 



In sections 4 and 5 we will look at the features of some of these solutions more closely. 



B. Class II Solutions 



We now consider fields in which 



II 



vld.l 



eS-^ = C + 



where A and C are constants. For A = we should get the interior Schwarzschild solution. 
Taking this relation and computing /(e'^) from ( p^ we have that 



and then (EG) tells us that 



fie'-) = eM~A^e-^n, 



-l^eM-A'e^n =F-DR'(l-^^^ 



1 

l4 



\n[-A'[F-DR'il-—) 



The electric charge, density and pressure of this solution turn out to be: 

Dr^ 1 



V2$ 



^Pi'') = ^ 



2^2 



D\ 



m 2$2 



[-ln(-A2$)] ^ 



-pM = -^ + ^[-M-A^$)] 



I^^'-2(l-^V$ 



Note that we must have < — A^$(r) < 1 to have real a solution. 



IV. PHYSICAL SIGNIFICANCE OF THE CONSTANTS 

In this section we relate the constants which appear in the previous solutions with the more physically relevant 
constants of total mass, total charge and radius of the body. To do this we consider the boundary conditions required 
of the solutions for them to match with the unique exterior solution. Before looking at some of these solutions in 
detail, we first draw a number of conclusions which come from their general form. 

Outside of the matter distribution we must have the Reissner- Nordstrom metric 



1 



2m 
r 



q 



If we let r = a be the boundary of the distribution, then mn, the total mass, and q, the total charge of the body are 
given by (Florides pq ) 



2 / y'^P 

'0 



' ' 2a 



Q{a) 



(22) 
(23) 



Our aim will be to relate the arbitrary constants which arose in the integration of the field equations to the physically 
significant constants m, q and a. We do this by using the continuity conditions which we expect the metric to exhibit 
across the surface r = a. 

In particular, we demand that e", e'^ and lu' all be continuous across r — a. (We note that as a consequence of 
the continuity of to' the pressure p vanishes at r = a by ( p^ ) - as expected). Hence 3 constants in each solution will 
be determined by the boundary conditions at r = a. Furthermore a 4th constant will also be expressible in terms of 
m, q and a by evaluating (p3) at r = a. Therefore, all of the constants appearing in the solutions presented (aside 
from C, the phase of the electric field, which we will ignore from here on) can be written in terms of the total mass, 
total charge and radius of the source. The actual relationships between these constants will obviously depend on the 
particular solution under consideration. However all the solutions have in common the Schwarzschild condition 



12 



which when put in (B3) gives 



Q2 3 

up -\ — = —r = constant^ 



i = ^ ("'-&■ <^^' 



This is equivalent to the continuity of e" across the boundary r = a. 

As a final comment we note that for the charge density to be non-singular at the origin, r = 0, we must have 

Q{r) ^ r'^ as r ^ 0, 

and this condition is satisfied by all of the new fields presented. 

A. Class I Solutions 

All of these solutions contain the constants D, F, R and either A or B. The constant R is given in terms of m, q 
and a by (p^, and A is positive for each of the solutions. 

We determine A and B in terms of the total mass, total charge and radius as follows. Recall that the charge inside 
a radius r is 



and given the functional relationship 



this becomes 



Q{r) ^ r^f e-i("+"\ 
Ae'^ =B + {C + ^f, 

Ar'^Lj'e-^^"-'^^ 
Q(r) = ± — . 

Now using (^3|) and the Reissner-Nordstrom solution we get that 

q 



which can be rearranged to read 



B^A 



'-|^)^a-^)(^f■ 



(25) 



For the Class la solutions B = Q and this gives A in terms of to, q and a. In this case we note that 

A>1 -^ m< \q\. 

Moreover, this relation means that 

A> 1 =^ -r < { — ] < ^■ 

A-\q) - 

Thus, A determines the balancing of the total mass and total charge. 
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For the Class lb solutions B ^ and A is 1 or i or |. Thus, equation ( p5| ) this gives B in terms of m, q and a. In 
this case we have that 

A>1 and B >0 => m < \q\, 

A <1 and B <0 =4> m > l^]. 

Turning now to the constant D, we note that putting our functional relationship in (|l8| ) gives 



/ i(w-a) 



W 6 2 



- = 2i:>Me'^ - B^. (26) 

r 

Since we require oj' to be continuous across r = a, we use the Reissner-Nordstrom solution and ( pq ) to get that 

a'' V a J \ a a^ 



1 / q 



2\ /„^ ..\ -2A 



Ai"^-^ 



(27) 

This will determine D in terms of the total mass, total charge and radius of the configuration for Class I solutions. 

The final constant F will be determined by the requirement that e" be continuous across r — a. This condition 
will differ for each solution and is rather complex, so we will take three illustrative examples in the Class I solutions. 

1. Class la A= \ 



Equation ( p5| ) yields 

'^1 -1 



2 / \ 2 

771 \ ^ 1 / "^ 9 



q J ^ \q a 



and so we have that m > \q\. We can solve this for m to get 



^2 



q , q 

m — - 
Equation (|27D reduces to 



777 = -^±^v/2(a2 + g2). 

a a 



a' 



The continuity of e'^ across r = a means F is given by 



/ 777 q \ I 2m Q \ 2 n , I 2777 o 

F=[--^\ 1 + \ +V21nl + \ 

q la J \ a a J \ a a'^ 



As a check, it's not hard to show that p(a) = 0, as required. 



2. Class la A = l 



We note that in this case the functional relationship ( |21| ) reduces to Majumdar's, and ( pSj ) tells us that 777 = \q\, as 
expected. The other constants work out to be 

a-i \ a J \q\ \ a 

To find F note that the pressure is 

^2^4 _ ^2 



p{r) = 



i?2$(r)2 ' 

and so for us to insist that p{a) = we must have F^ — D^R^. Thus the matter pressure vanishes, and we are led to 
the charged dust solution first discovered by Cooperstock and De La Cruz [El| . As this solution has been considered 
before, we will not discuss this solution further here. 
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Class lb A^l B <0 



The constants B, R and D are given by 



B = l-(™ 



1 q /2m q 

R^ a^ \ q a 



w^ \ a 



Obviously 



m > \q\ <=^ B <0. 



We shall restrict our discussion to the undercharged case B < (i.e. m > |g|). In this case the remaining constant F 
is given by the continuity of w' across r = a: 



sec 



- 1 \ F- DR^ 1 - 



-.2 \ 2 



i?2 



m q 
q a 



which we can rearrange to read 



F 



q 



m q 
q a 



q 



1- 



2m q^\^fm q\ f 2m q^ 



a a^ 



q a 



q a 



(28) 



B. Class II Solutions 

These solutions contain the constants A, D, F and R. Equation (|2J) gives R in terms of m, q and a. 
We determine A in terms of the total mass, total charge and radius as follows. Recall again that the charge inside 
a radius r is 



= r-2 ,/ -i(a+^) 



Q{r) ~ r ip'e 



and given the functional relationship 



this becomes 



:^Aei-^C 



Q{r) 



Ar^ 

2V2 



tj'e-2"+4'^. 



Now using (E3|) and the Reissner-Nordstrom solution at r = a we get that 



A=V2i'^-^ 
q a 



2m q^ 



Turning now to the constant Z?, we note that putting our functional relationship in (Rq) gives 
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^/gi(w-a) 



2L»exp 



(A^e^-) 



Since we require uj' to be continuous across r = a, we again use the Reissner-Nordstrom solution to get that 



D = -3 ( m - ^ ) exp 



-21™-^ 



1- 



2to q^ 



(29) 



(30) 



q a ^ 
The final constant F will be determined by the requirement that e" be continuous across r ^ a. This works out to 



be 



F 



q a 



. exp 



-2(^-^ 
q a 



m q 
q 2a 



cj(a) 



"1 q \ ^-u(a) 

q a 



where e'^^"^ = 1 - 2m _^ 9^ 



V. REGULARITY AND ENERGY CONDITIONS 



In order for a solution of the field equations to be considered physically plausible there are a number of conditions 
that it must satisfy. In particular, we may require that a solution: 

(a) be non-singular, 

(b) have no horizons, 

(c) satisfy regularity conditions at the origin, 

(d) satisfy energy conditions. 

The first of these simply means that the metric and the energy- momentum tensor are finite throughout. This is 
not true for the charged dust solution (Class la, A= 1) first discovered by Cooperstock and De La Cruz |£l|, and we 
exclude this solution from our consideration. We will consider the other solutions shortly. 

A horizon will exist where e'^ = 0. This can happen either outside or inside the matter distribution. Outside the 
distribution, horizons can only occur if ?7i > \q\ and then these can occur at 



r± = m zL \/vfi^ — q2 



However, from equation (p3) we see immediately that for all of the solutions 



m > — . 
2a 



We also note that our derivation of the solutions has implicitly assumed that lo' > {cj. equation (18)). Since we 
require that uj' be continuous across r = a we must have 



m > 



(31) 



We therefore conclude that there will be no Reissner-Nordstrom repulsion outside of any of the solutions which we 
have derived. (For a full description of this phenomenon see, for example, Papapetrou [p7|). 
Furthermore, it is not hard to show that 



< 



r 

m 



< 



'+j 



and since we cannot have the boundary of the body lying in the non-static region between the horizons, we conclude 
that for all of the solutions presented, if m > \q\ holds, then (by equation (pl|)) 
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a > m + \/m?- — q^ , 

i.e. there are no horizons exterior to any of the solutions. As for the possibihty of horizons inside he matter, we will 
discuss this shortly. 

The requirements that the solution be regular at the origin are 

(i) a'(0)=0, 

(ii) u'iO) = 0, 

(iii) p'(0) = 0. 

A straightforward differentiation of equation (O) shows that (i) holds for all interior solutions satisfying the 



Schwarzschild condition (|16|), while equations ( |26|) and (29) ensure that regularity condition (ii) holds for all Class I 



and Class II solutions which are non-singular at the origin. Below we will consider condition (iii) for the new solutions 
individually. 

Finally there are many different energy conditions that can be imposed on the solutions. We consider: 

i) -Ti > 0, 

ii) Tl — T^ > or p + p > the weak energy condition, 

iii) T" — T4 > the strong energy condition. 

We could add to these the extra conditions 

iv) p > 0, 

v) p > 0, 

vi) 3p + p > 0, 

vii) p' > 0, 

viii) p' < 0. 

Of these, condition (i) is automatically satisfied by all solutions (c/. equation ([16|)). Note also that energy condition 
vi) implies the strong energy condition (iii). To continue our investigation we must look at the explicit solutions 
separately. We will see that all of the new solutions presented satisfy (i) to (iii), while many of the conditions (iv) to 
viii) also hold in some cases. 

A. Class la A = i 

This solution is clearly non-singular everywhere and has no horizons. 
Also, since 

kT^ - kT* = D'^r* -I- 3\/2 D ( 1 - -^ y > 0, 

1 

3kp + kp= ii:>V^ + 3\/27:> f 1- ^ ) >0, 



i?2 



energy conditions (ii), (iii) and (vi) all hold. 

Differentiating the expression for the matter pressure gives 



Kp' — Dr 



\/2 / r^\ ^ 



*^-i^l^-i?^j 
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Thus, p'{0) — and the solution satisfies all of the regularity conditions at r = 0. Focusing on the term in square 
brackets, and using the expressions for D and R in terms of rn, q and a, as well as the fact that m > \q\ 



hD 



^(i 



i?2 ^ i?2 



<iD 



V2 



2a 

^2 



m 2\f2 I q 

< — :r 7T- \ fn — — 

2a 



'2\/2-l' 



V2 q^ 
2%/2- 1 a 



<0, 



where this last inequality comes from equation (|3l|). Thus we conclude that p' <0 (energy condition (vii)) and since 
p(a) = 0, p > (energy condition (v)). 
Finally, differentiating the matter density 



Kp' = -^D^r < 0, 



satisfying energy condition (viii), and since 



Kp{a] 



3 D^a^ 
~ i?2 ~ 4 

6 
a 
>0 



2q2 



3 1™ Q 

■^ \ da 



(32) 



we have p > (condition (iv)). 

In conclusion, this solution satisfies all of the regularity conditions and energy conditions (i) to (viii) . 



B. Class lb ^ = 1 B < 



For this solution to be non-singular and have no horizons, we require that 



< 



1^1 -1 

q 



$(r) < |. 



Since $(r) is a non-decreasing function of r, this is equivalent to 



$(0) > and $(a) < -. 



1^1 -1 



The second of these follows from the definition of the constant F, equation (p8|). The first can be seen to hold as 
follows. 

A short calculation shows that 



2m q^ 

1 + -^ >0 ^ 

a a^ 



^1 -1 



> 



m \q\ 



(33) 



Now, 
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$(0) ^ F - DR^ 

2 

m 

q 



m q 
q a 



+e^ 



la;(a) /^™ _ q\ [^ _ g 



> 



> 



- 1 



rn q 



m q \ I Zm q 



2m q 



m \q\ 



q a J \ q a 
TO_|g|\" 1 2711 \q\\^ 



>0, 



where we have used equation ( |33D on Hne four. Thus the solution is non-singular and has no horizons. 
A straightforward calculation shows that p'(0) = 0, and therefore the solution is regular at the origin. 
The charge density and electrostatic potential turn out to be 



KCT = ±2D 



1 $(r) 





{<-i 


1 

) +Dr^ 
J 




_\qJ 


1 
2 

sec 


7/ 



^1 -1 

q 



tan 



-) -1) $(,) 



-) -1) ^r) 



Thus, the charge density does not change sign throughout the distribution of matter, and 



3kp +Kp = 2D^r^ 



1^1 -1 

q 



sec 



-) -1) *w 



A2D 



1 $(r) 



2 \ 2 



i?2 ) 



± 



e 2-Ka 



des" 



-AvCr, 



in agreement with equation (H). 

Finally, turning to the energy conditions, we can see from the above calculation that 3p + p > and so the strong 
energy condition holds, while 



Kp + up = — - + 2y/—BD sec 
R^ 



-B$(r) 



-B$(r) 






so the weak energy condition also holds. 
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C. Class II 



For this solution to be real and have no horizons we must have 

_i_ < $(,) < 0. 

Since $ is an increasing function of r , this will be true as long as 



*(a) < 



and 



$(0) > 



^2- 



The first of these is trivially true, while the second is equivalent to 



1 / \ -1 

m q \ m q 



2 2 

m \ m q 



q a J \ q 2a 
A short calculation shows that this holds if 



m _^ ^_ i.(a) 



exp 



q a 



> -1. 



-I I Hi 1 in 



2a ' 2a2 



,^uj(a) 



^ I g i 2a ^^ 2q2 <i ' ' 



exp 



q a 



< 1. 



Since both of the factors on the left-hand side are less than 1, this inequality holds. Thus the solution has no horizons. 
Differentiating the pressure, it can be shown that p'(0) = and the solution is therefore regular at the origin. 
Finally 



-^ + -p = f-J^[-inM^*W)]-^(i-^ 



2 \ 2 



>0 



3kp + Kp = 



D 



$(r)2 
>0. 



[-\n{-AH{r))] ^ 



Dr' - 6$(r) ( 1 - ;^ 



We conclude that both the weak and strong energy condition (as well as condition (vi)) hold for this solution. 

VI. CONCLUSIONS 

Functional relationships between the gravitational and electric potentials can take many forms. In the presence 
of matter pressure, these relationships can go beyond the original quadratic form introduced by Weyl, and can give 
rise to matter distributions that satisfy many physically reasonable criteria. In addition these fields satisfy a relation 
between their matter, electric field and charge densities. This generalises the fact that the active gravitational matter 
density and charge density of a Majumdar-type field are numerically equal. For Majumdar-type fields the total mass 
and charge as measured at infinity are equal, while bounded non-singular perfect fluids with non-vanishing pressure 
cannot generate such fields. 

The author would like to thank Petros Florides, under whose supervision most of the above work was carried out, 
and the referees for many helpful comments and suggestions. 
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